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Abstract 



A relativistic constituent quark model is adopted to give an unified description of the leptonic and 
semileptonic decays of pseudoscalar mesons (vr, K, D, Da, B, Bs). The calculated leptonic decay 
constants and form factors are found to be in good agreement with available experimental data and 
the results of other approaches. Eventually, the model is found to reproduce the scaling behaviours 
^vj \ of spin-fiavor symmetry in the heavy-quark limit. 



> ; 1 Introduction 

Semileptonic decays of pseudoscalar mesons allow to evaluate the elements of the Cabibbo-Kobayashi- 

Maskawa (CKM) matrix, which are fundamental parameters of the Standard Model. The decay K -^ nev 

provides the most accurate determination of Vus, the semileptonic decays of D and B mesons, D — > 

Q^ . K{K*)lv, B —t D{D*)lv and B -^ ■n{p)lv, can be used to determine |Ks|, \Vcb\ and |14&|, respectively. 

0^ ' This program can be performed if the non-perturbative QCD effects, which are parameterized by the 

ci ■ form factors, are known. Up to now, these form factors cannot be evaluated from first principles, thus 

Q^l models, more or less connected with QCD, are usually considered for this purpose. Here we discuss a 

relativistic quark model pi , previously used to study the baryon form factors B . 
flj " This model is based on an effective Lagrangian describing the coupling of mesons with their constituent 

r^ , quarks. The physical processes are described by the one-loop quark diagrams and meson-quark vertices 

related to the Bethe-Salpeter amplitudes. In principle, the vertex functions and quark propagators should 
be given by the Bethe-Salpeter and Dyson-Schwinger equations, respectively. This kind of analysis is 
provided by the Dyson-Schwinger Equation (DSE) studies Q and an unified description of light and 
heavy meson observables was carried out in fj, |5[|. Here, instead, we use free propagators for constituent 
quarks and consider a Gaussian vertex function as Bethe-Salpeter confining function. The adjustable 
parameters, the widths of Bethe-Salpeter amplitudes in momentum space, and the constituent quark 
masses, are determined from the best fit of available experimental data and some lattice simulations. 
Our results are in good agreement with experimental data and other approaches. We also reproduce the 
spin-flavor symmetry relations and scaling for leptonic decay constants and semileptonic form factors in 
the heavy-quark limit [pi . 

An approach similar to the one presented here, can be found in M. That model is based on meson- 
quark interactions, so the mesonic transition amplitudes are described by diagrams with heavy mesons 
attached to quark loops. The free propagator has been used for light quarks, while the quark propagator 
obtained in the heavy-quark limit has been adopted for heavy quarks. 

2 Our model 

Our starting point is the effective Lagrangian describing the coupling between hadrons and quarks. The 
£int{x) ^ gHH{x) dxi dx2^H{x;xi,X2)q{xi)rHXHq{x2) (1) 

^Based on talk given by the first author at the QFTHEP '99 Workshop, Moscow, May 27- June 2, 1999. 



describes the transition of the meson H (Xh is the corresponding combination of Gell-Mann matrices, 
Th are Dirac matrices) into its constituents qi and (72- The function ^h is related to the scalar part of 
the Bethe-Salpeter amplitude. 

The coupling constant gn is given by the derivative of the meson mass operator Hh by the compos- 
iteness condition H: 



Zh = i- ^n'^{<) 



0. 



(2) 



It is worth noticing that, due to the absence of confinement, the sum of constituent quark masses 
should be larger than the mass of the corresponding meson otherwise, imaginary parts in physical quan- 
tities appear. This allows us to consider low-lying pseudoscalar mesons only. 

Now, to give an example of the hadronic part of invariant amplitudes, we will evaluate the form factor 
/+, which appears in the semileptonic decays of a pseudoscalar meson into another one, H — > H'iv. The 
invariant amplitude can be written as: 



A{H{p)~^ H'(j>')ev) = 



^I^[i^^il-^,)^]MU,{p,p'l 



(3) 



where Gp is the Fermi constant, V^^g^ is the corresponding Cabibbo-Kobayashi-Maskawa matrix element 
and, in our model, Mf^^,{p,p') is given by: 



4 7r2 ./ 47r2^ 

igngH 



,(j>H{-k')^H'{~k')tr 



rSs{^h'S2{^+ /)r{l-To)Si{^- 



(4) 



4 7r2 



U{p^p",q'){p + pr + Up',p'^q')q' 



Here, q — p — p' , Si{k) ~ l/{mi— K) is the propagator of the quark i with mass nii 



To evaluate the integral in Eq. (|4|) we have to calculate the following integrals {J-{—k^) = (f)H{—k^ 

^H'i-k')): 

cf^k {l,k^',k^'k'',kl'k''k^)T{-k'^) 



j{0,fi,i.ii^.fiL'5) __ 



(5) 



TT^i [ml — {k+ p)2] [m2 ~ {k + p')^] [m| — k'^] 
They can be evaluated using the standard Feynman a-representation and the integral Cauchy represen- 
tation for the function J^: 



, . j.i^j_jj,.^i^.^t.:)^yiM^ 



(6) 



J^ = - dt 



l + t 



d^a5[ 1 -^a, 



pM 



dT{zi) 
dzi 



(7) 



n * — -L 



(8) 



n 1 — -L 



M^p5 _^ nSpU _^ ^5pM 



l+t 



rn^i) 



piJ. pv pd 



l + t 



dzi 



(9) 



where P = aip + a2p' , -D3 = aia^p^ + a20Lzv''^ + aia29^, and 

3 



Zl = t 



y^ oiim\ - Da 



,1=1 



(1 + 



p. 



(10) 



Finally, we have the analytical expression for /_| 



/+(p^p'^?') 



4-3(ai+a2) 



+ 1 -^ I (mi + m2)m3 + I -(ai + a2)[mim3 + 17121713 — mim2J 



+aip^ + a2V - P 



l + t 



2 - (ai +a2) 



1 + i 



(11) 



which can be used also to obtain the normalization constants gn and gn' in Eq. (g). For example, gn is 
given by 



47r2 



5H 



3/+(p2,p^0) 



(12) 



where we put mi — 7712 = m. 



3 Results and discussion 

Once obtained the analytical expressions for the invariant amplitudes, a comparison with the experimental 
data is in order. For this purpose we have to specify the analytical form of (/>// (— A:^) and the constituent 
masses appearing in the expressions. In particular, we choose the Gaussian form (f){—k'^) — exp{fc^/A|^} 
in Minkowski space, where the magnitude of Ah characterizes the size of the BS-amplitude and is an 
adjustable parameter in our approach. Thus, to describe processes involving tt, K, D, Dg, B, and Bg 
mesons, we have six A-parameters plus four quark masses (tti^ — 111^ — mj_, m^, 771^, and mi,), all of 
them are fixed via the least-squares fit to the observables measured experimentally or taken from lattice 
simulations (as is reported in Table |l|). 



Table 1: Calculated values of a range of observables {gyrj-y in GeV~^, leptonic decay constants in GeV, 
form factors and ratios are dimensionless). The "Obs." are extracted from Refs. [g|, |l^. 111, 
The quantities used in fitting free parameters are marked by "*" . 
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IB II- 







Obs. 


Calc. 






Obs. 




Calc. 


* 


9-1^71 


0.274 


0.242 


/f-(0) 




0.98 




0.98 


* 


u 


0.131 


0.131 


* /f^(0) 




0.74 ± 0.03 




0.74 


* 


Ik 


0.160 


0.160 


/r(o) 








0.73 


* 


Id 


0.19lt^^ 


0.191 


/f-(0) 




0.27 ± 0.11 




0.51 


* 


Id, 
fo 


1.08(8) 


1.08 


Bi{K - 


■nil') 


(4.82 ±0.06) • 


10-2 


4.4- 10-2 




Id. 


0.206l^^ 


0.206 


Bt{D -^ 


Klu) 


(6.8 ±0.8) -10 


-2 


8.1-10-2 


* 


fs 


0.1721^7 


0.172 


Br(B^ 


Dlu) 


(2.00 ±0.25) • 


10-2 


2.3-10-2 


* 


fB., 
fB 

Ib. 


1.14(8) 


1.14 
0.196 


Br(B^ 


nlv) 


(1.8 ±0.6) -10 


-4 


2.1- 10-4 



The fitted values for the free parameter of our model are the following: 

(A^, Ak, A^,, Ad,, As, AbJ = (1.16, 1.82, 1.87, 1.95, 2.16, 2.27) GeT^ 



(13) 



(m,, TOs, TOc, rrib) 



(0.235, 0.333, 1.67, 5.06) GeV . 



(14) 



Note that the A values are larger for mesons with larger mass, i.e. Ah < Ah' when mn < mn'. 
This correctly corresponds to the ordering law for the sizes of bound states. 

The u(d)-quark mass and the parameter A^ are almost fixed from the rate tt — > fiv and tt'^ —> 77 with 
an accuracy of a few percent. Moreover, the obtained value of tUu is less than the constituent-light-quark 
mass used in baryon physics. 

Let us now consider the g^-behaviour of the form factors. Since a numerical integration should be 
done (see, Eq. (O)), we do not have a simple analytical expression for them. However, looking at the 
numerical results, the form 



f^ il') 



1 



ml 



q 



'H 



(15) 



is suitable for a good description of the g^-behaviour, once the parameters 60, &i and /(O) are fixed. 
Their values are collected in the following Table: 



(16) 



It should be noted that a value for /:^'^(0) larger than those_obtained by QCD Sum Rules |16|, and 





K ^TT 


D -^ K 


B ^ D 


B-f-K 


/(O) 


0.98 


0.64 


0.73 


0.51 


bo 


0.28 


0.64 


0.77 


0.52 


hi 


0.057 


0.20 


0.19 


0.38 



other approaches |18 . In any case, as one can see from the Table ^ the agreement between our predictions 
and experimental data on semileptonic branching ratios is impressive. 

As we have seen, our model gives an accurate and unified description of the weak and radiative 
(tt" -^ 77) transitions involving pseudoscalar mesons. Moreover, as already stated in the introduction, it 
is able to reproduce the scaling behavior predicted by QCD in the heavy-quark limit. For more details, 
see the original paper. Here we report the way to obtain this limit in the expression for /_|_. The heavy- 
quark limit corresponds to consider mi = M — *■ 00, TO2 = M' — > 00 and p^ = {M + EY , p'"^ = (M' -I- EY 
with E being a constant value independent of M and M' . By replacing in Eq. (|l^) the variables ai with 
ai/M and a2 with a2/A/', one obtains 



where 



z = u-2E 



M + M' 
2MM 







^^.Jdtl^j^j Jdaajdr(^-T'{z) 



00 

/ du!F{z 



1 00 
M + M' If dr 



2MM' 2 J W 




v/^W^ 



m'' + z 



w 



W = 1 + 2t(1-t)(w-1), 



at 

YTt 



2MM' 



(17) 



(18) 



Therefore, using the relation between /_)_ and /+ and the normalization condition, the correct scaling 
relation is found: 



f+ -^ n'.^"} c(w) oc / — / du(hf,(z] ^ 

■'^ 2VMM^ ^ ' J Wj ^"^ ' m^ + z 





(19) 



In conclusion, we can see that the agreement with experimental data and lattice results is very good, 
with the exception of the value of /^" (0) which is found to be larger than the monopole extrapolation of 
a lattice simulation, QCD Sum Rules (cf. |]16|) and some other quark models (see, for example, |l7[ [iq]). 
However, this result is consistent with the value calculated from Refs. g, |l^ and a light-front constituent 
quark model [E0|. Moreover, it allows us to reproduce the experimental data on i? ^ nlv decay. 
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